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Abstract 

This note is concerned with an important for modelhng question of existence of solutions 
of stochastic partial differential equations as proper stochastic processes, rather than 
processes in the generalized sense. We consider a first order stochastic partial differential 
equations of the form ^ = DW, and ^ ~ ^ = DW, where Z) is a differential operator 
and W{t, x) is a continuous but non-differentiable function (field). We give a necessary and 
sufficient condition for stochastic equations to have solutions as functions. The result is 
then applied to the equation for a yield curve. Proofs are based on probabihty arguments. 

1 Introduction 

Stochastic differential equations are often obtained from ordinary ones by introduction 
of noise, which is taken to be a generalized function (Schwartz distribution). Following 
Rozanov 1995, the noise in SPDE is defined as follows. Let W{t,x) be a continuous 
function (field) and D a differential operator. If W is not differentiable, then DW does 
not exist in the usual sense of a function, but can be considered as a generalized function. 
It is defined by the following action on any test function (p{t,x) G C^((0, oo) x (0, oo)), 
the space of infinitely differentiable compactly supported functions, 

{DW,ip) = {W,D*ip) = J W{t,x)D*ip{t,x)dtdx, 

where D* is the adjoint operator of D. The adjoint operator is defined by the identity 
(Df, (p) = (/, D*(p) that holds for all smooth functions / and test functions (p. For details 
see Rozanov 1995 p. 99-103, or Lang 1993. 
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In applications the noise is typically random. White noise, for example, is a generalized 
derivative of Brownian motion B{t) = and the noise considered in second order 

partial differential equations is the second derivative of Brownian sheet B{t, x), 13{t, x) — 
^{t,x), (Walsh 1984, Carmona and Nualart 1988, Freidlin 1988.) First we recall a 
definition of Brownian sheet B{t,x) on positive quadrant of the plane R"^, which is the 
basic model for White noise field with its various modifications. In our applications t is 
time and is non-negative, as well as the state variable x. This is not really a restriction, 
if needed the Brownian sheet can be taken on the whole plane i?^. A Gaussian random 
measure B on R"^ is defined by the following properties. For any Borel set A C R\, B{A) 
is a Normal random variable with zero mean and variance given by the area of A. For non 
overlapping Ai, A2, B{Ai) and B{A2) are independent and B{Ai U A2) = B{Ai) + B{A2). 
Put B{t,x) = B{[0,t] X ([0,a;]). It is known that for almost all realizations, B{t,x) is a 
continuous but nowhere differentiable function of t, x. For details see e.g. Walsh 1984. 

In the next section we give results for solutions of first order partial differential equa- 
tions to be functions when the noise is a generalized function, and then use probabilistic 
arguments to derive conditions for the case when noise is obtained from a stochastic pro- 
cess, such as a white noise. In Section 3 we apply the result to the equation of a yield 
curve. Section 4 contains the proofs. Since our focus is on modelling with SPDE's it is of 
prime concern that solutions will be proper functions. The classical existence and unique- 
ness results apply in the space of generalized functions. Perhaps a classical approach is to 
use existence and uniqueness in the space of generalized functions and then apply some 
regularization results. Our approach is more direct, and utilizes the probabilistic nature 
of noise in the equations. 



2 First Order SPDE's 

Let D be a first order differential operator 

Dip{t, x) = a{t, x)^it, x) + b{t, x)^it, x) + c{t, x)f{t, x), (1) 

where a and b are smooth functions from C^, and c is continuous. The adjoint operator 
D* is easily found to be 

at ox 

It turns out that only a particular form of the differential operator D yields solutions as 
functions of the SPDE (2). 

Theorem 1 Let D is given by (1) and U be a solution to the pde 

f = (2) 
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in the sense that — (U, = {W, D*ip) . for any test function cp. Then it holds 
[Uit,y)-U{Q,y)]dy 

= /* [bis, x)W{s, x) - b{s, 0)W{s, 0)] ds + r [a{t, y)W{t, y) - a(0, y)W{Q, y)] dy 
Jo Jo 



+ 



f I 

Jo Jo 



dy. 



(3) 



Thus a solution of (2) as a function exists if and only if / b(s,x)W(s,x)ds is differ- 

Jo 

entiable in x. In particular, when b is identically zero, a solution as a function exists. 



Consider now the equation 



dr dr 



= DW. 



dt dx 

We say that r{t, x) is a solution of (4) if (5) holds for any test function (p 



(4) 



J r{t,x) 



'dip d(p^ 



{t,x)dtdx = J W{t,x) 



f d{a(p) d{b(p) 



dt dx ) ^'-^ J' ■ — V dt + ^ - '^J 



Theorem 2 Suppose that the noise process is such that for any x > 0, W{0, x) — 0. 
Then the equation 

dr dr ( dW dW \ 

-g^{t,x) - —it,x) = a{t,x) {-g^it,x) - -g^it,x)\ + c{t,x)W{t,x) (6) 

(this is equation (4) with b{t,x) = —a(t,x)) with the initial condition r(0,,x) = ro{x) and 
r{t, x) is continuous at t = 0, for all x, has for unique solution the function 



r(t,x) = a(t,x)W(t,x) + W(s,x) 

Jo 



da da 

^(S'^J - -Q^is,x) + c{s,x) 



ds + ro{t + x). (7) 



Next (motivated by applications) we take the noise W{t,x) as the Brownian sheet 

W{t,x) ^ B(t,t + x) -.^ B'^it). (8) 
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Theorem 3 A solution to equation (4), where D is given by (1) and W{t,x) by (8), as 
a function exists if and only if 

a{t,x) ^ -b{t,x). (9) 
Moreover, in this case, the general solution ( 7) can be written in the form 



r{t,x)^ f\{s,x)dB''{s)+ t B'^is) 
Jo Jo 



da 
dx 



{s,x) + c{s,x) 



ds + ro{t + x), 



where the integral Jq a{s, x)dB^{s) is the usual Wiener-Ito integral. 



3 Application to Yield Curves 

Let r{0,x) = ro{x) be a given yield curve at time t — 0, that is the interest on the 
investment maturing at x. Let r{t, x), t,x > denote similar curve at time t, that is the 
interest rate at time t on the investment maturing ai t + x. In the absence of noise it 
should hold that r(t,x) = r{0,t + x) = ro(t + x), otherwise one can make a riskless profit. 
Assuming that ro(a;) is smooth, the evolution of the yield curve therefore is described by 
pde 

dr dr , ^ 

with the initial condition r(0,a:) = ro(a;), see Musiela and Sondermann 1994. Consider 
now a stochastic analogue of this pde given by (4) || — |^ = DW., where x) = 
B{t,t-\- x). The reason we take B{t,t + x) rather than B{t, x) is because we model the 
yield at the point t+x, and it would be natural to write r{t, t+x) instead of r{t, x) and add 
to it the basic noise at the same point {t,t + x). The process W{t, x) — B{t,t + x) = B^{t) 
for any fixed a; > 0, is not a Brownian motion, however it is a continuous martingale with 
independent but non-stationary increments. Condition 1^(0, a;) = for any a; > in 
Theorem 2 means here that at time zero there is no uncertainty in the yield curve. 

Musiela and Sondermann 1994 introduced the noise by considering for all x > the 
stochastic differential equations 

dr{t, x) = a{t, x)dt + a{t, x)dW{t), (11) 

which hold for all a; > 0. In comparison with this equation, the general solution (7) for 
any fixed x > satisfies 



dr{t,x) = a{t,x)dB''{t) + 



B'^it) (^{t, x) + c{t, x)) + r'oit + x) 



dt. 



Hence in the spde solution the "drift" term includes random noise, moreover the equations 
are driven by martingales B^{t) dependent on x. 
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4 Proofs 



The proof of Theorem 1 requires the following basic lemma. 

Lemma 4 Let g be an integrable function on [0, +oo) x [0, +oo), and assume that for any 
test function, on (0, +oo) x (0, +oo), (f, 



5(^,2;)^^(t,x)cix = 0. 
Then, g{t, x) = g{t, 0) + ^'(0, x) + const. 
Proof of Theorem 1 

For any test function (p, — (U,^^ =< W, D*ip >, or equivalently 

d{aip) d{bip) 



(12) 



dt 



dx 



— C(p 



(t, x)dtdx 



I U{t,x)^{t,x)dtdx = J W{t,x) 

= fw{t,x) 
Integrating by parts, we get 

rU{t,y)dy- r a{t,y)W{t,y)dy- f' b{s,x)W{s,x)ds+ 
Jo Jo 



d(p ^dcp f da db \ 
dt dx \ dt dx ^ I ^ 



{t, x)dtdx. 



(3) now follows from Lemma 4. 



dy 



d\ 
dtdx 



{t, x)dtdx ~ 0. 



□ 



Proof of Theorem 2 

Make a change of variables t — t, ^ — t -\- x, p{t, ^) — r{t, x), Q;(r, ^) = a{t, x), P{t, ^) = 
b{t,x), 7(r, ^) = c(t, x), iI^{t,^) = (p{t,x) and V{t,^) = W{t,x). Perform the above 
change of variables in (5) and use the relation b{t,x) — —a{t,x) to obtain 



(13) 



It follows from Theorem 1 that 



« r (da 



[p(T,77)-p(0,77)]d77 = air,rj)V{r,rj)dri + | t::---c]{e,ri)V{e,rj)dedri. 
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Differentiating with respect to ^, we get 

p{r, = OV{r, - [ V{9, - 7) {0, Ode + p(0, 0- (14) 

Going back to the original variables and taking into account the initial condition we obtain 
(7). Notice that a formal derivation of (14) is obtained by integrating with respect to r 
the equality ^ — + (to be understood in the sense of generalized functions). □ 



Proof of Theorem 3 

The proof of Theorem 3 is based on the following lemmas. 

Lemma 5 Let X be a Gaussian measure on a measurable space {E,S) with intensity [x, 
R and S be two functions in L^{E, £, /i), F and G be two measurable sets, and {F^)k and 
{G^)k be measurable partitions of F and G respectively. Assume that 

1. fJ>{F) < +00 and fJ,{G) < +00; 

2. lim„^<^supi<fe<„/x(FjJ^) = lim„^<^supi<fc<„/x(G^) =0. 
// {F^ n Gl)k is a partition ofFnG, then 

n „ 

in the L'^-sense, where {R^) and {S^) are approximating sequences on {F^ n G1)k of R 

and S respectively. 

If, on the other hand, for all k and I, one or both of F^ n and Fp n G^ are empty (in 
particular F^ nGl = i), then 



y^^Y.Rlsix{F^)x{Gi)^Q 

in the L'^ -sense. 



n— >oo ■ 



Proof of Lemma 5: 

Let us first note that Yax{X{Fl')X{Gl)) = ii{Fp)^i{Gl)+ ii{FlpnGlf , and that, for A; ^ 

coy{x{fj:)x{GI),x{f-)x{G^)) = ^i{F]:r^G^)^i{F-r^Gl). it follows that, if {Fj:nGi)k 
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is a partition of F n G, then the random variables X{F^)X{Gl) and X{Fp)X{Gf), for 
k I, are independent. 



k=l -^^nG 



< 



k=l k=l 



+ 



jk=l 



FnG 



RSdfi 



The second term of the right hand side goes to zero. By the independence of the random 
variables X{F^)X{Gl),th.e first term becomes 



ye.r[Y,R]:S-,X{F;:)X{G-,) 

\k=l J 



= EiRkSkrye.r{X{F-)A:{G^^)) 

k=l 

n 

= j:iRkSkfHKMGl) + fi{F;:nG: 

k=l 

n 

< ^T.iRkS^fi'iKuGiy 



n\2i 



fe=l 



< 2 sup ^l{F^UGl)Y.{RlS-,f^^{F;:uGl), 



l<k<n 



k=l 



which completes the proof of the first statement. The second statement follows in the 
same way. □ 

Lemma 6 Let X be a Gaussian measure on a measurable space {E, £) with intensity ii, 
F be a measurable set, and (-F^)fc be a measurable partition of F. If there exists k > 
such that 

lim n'^ sup uiF^) = 0, 



then, with probability one. 



lim sup \X{F-)\=0. 



Proof of Lemma 6: Fix £ > 0. 



sup \x{f;:)\>s 

l<k<n 



1-P 



sup |A'(F«)|<£ 

l<fe<n 



< 1-n 

fe=i 



< 1 



1 — exp 
1 — exp 



where ^{F"-) = supj^<j^<„ The proof is ended by the fact that 

^ s 



Eh 



1 — exp 



< +00. 



□ 



Proof of Theorem 3: The " if " part is given in Theorem 2. We now prove the " only 
if " part. Using the notations of the proof of Theorem 2, we obtain, instead of (13) and 
without the assumption b{t,x) — —a{t,x), 



dtp 



difj f da d{a + (3) 



-7 U; 



It now follows from Theorem 1 that for a solution to exist in the space of functions, 

Zit.x) := f Ais,x)Bis,s^x)ds^ A(s,x)W(s,x)ds ^ r (a + (3)(e,C)V(e,C)d0, 
Jo Jo Jo 

where A = a + b, must be differentiable in x. Let Y{t,x) = A{t,x)B{t,t + x), x < y, 
S'^ — [y — x)/n and 5^ = M". We have, 

j:(z{t,x + 5l)-Z{t,x + 5U)" 

k=l 

= E / {Y{s,x + 6l)-Y{s,x + S]:_,))ds 

ft ft ^ 

= / j:[{Yi^^^ + ^k)-Y{r,x + S]:_,)){Y{s,x + Sl)-Y{s,x + S]:_,)) 

•"^ -^0 fc=l 



drds. 



Now AY,-{s) := Y{s,x + 6-,)-Y{s,x + 6^^^_,) = Al{s)B{F^{s)) + BU{s)A{Gl{s)), with 
Al{s) = A{s,x + 5l), Bl_,{s) = 5(5,5 + a; + 5^), F]:{s) = [0,5] x[s + x + 5l_,,s + x + 5ll 
G^{s) = [0, 5] X [x + x + 5^], S is the Gaussian measure associated to the Brownian 
sheet B, and A the measure associated to the function A. Thus 



EAn"(r)An"(5) 



k=l 
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E Kir)Al{s)B{F-{r))B{Fj:{s)) + Al{r)BU{s)B{F^{r))A{Gl{s)) 



k=l k=l 

n 



k=l k=l 

Now each one of the last three terms converges, with probabihty one, to 0. Indeed 



k=l 



< sup \B{F^{r))\ 

l<k<n 



Al{r)BU{s)A{Gl{s)) 



k=l 



which, according to Lemma 6, goes to 0. Recall that, if /i is the intensity of the Gaussian 
measure B, then 

sup ^i{F^{r)) = ry^. 

l<k<n n 

The same goes for the third term. The convergence to of the fourth term follows from 

lim sup \A{F^{r))\ ^ 0. 

l<k<n 

To investigate the asymptotic behavior of the first term, we first make the following 
observation. For r ^ s, there is no {no > {y — x)/\r — s\) such that for any n > no, and 

any k and /, FJ}{r) fl Fp{s) = or Fp(r) fl FJ^^s) = 0. Applying Lemma 5, we see that 
'Ek=i^Ur)Alis)B{F^{r))B{F^{s)) converges in to 0. On the other hand, if r = s, 
then, by Lemma 5, 



jim ± Al{sfB{F]:{s)f = r A{s, zfsdz 



k=l 

in the sense. It follows that 



n— >oo 

fc=l 



in probability, and that 



i™o E {Z{t, X + 51) - Z{t, X + 5U)" ^ fir A{s, zfsdz) ds (15) 

k=\ \'' X / 

in probability. Now assume that Z{t,x), as a process in x, is locally Holder of order 
£ > 1/2, then 

Z{t, a; + 5^) - Z{t, x + | < const 

and 

n 

{Z{t, X + kSn) - Z{t, x + {k-l)Sn)f < const^ {y - x) — ^ 0, as n ^ oo. 

k=l 

This combined with (15) shows that A{s,x) must be nil. □ 
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